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ABSTRACT 

£^ . We study zero temperature correlation functions of the spin-^ Heisenberg XXZ 

model in the critical regime — l<A<lina magnetic field by means of the 
Dual Field Approach. We show for one particular example how to derive determi- 
nant representations for correlation functions and how to use these to embed the 
correlation functions in integrable systems of integro-difference equations (IDE). 
These IDE are associated with a Riemann-Hilbert problem. 



1. Introduction 



The evaluation of correlation functions in integrable one-dimensional quantum 
systems is one of the main outstanding problems in Mathematical Physics. Quite re- 
^ ■ cently there has been exciting progress in this direction: the Kyoto group succeeded 

in deriving integral representations for some correlation functions of the Heisenberg 
XXZ model in the massive regime A > 1 by takine advantage of the infinite quantum 
affine symmetry of the model on the infinite chairfM. These integral representations 
are most powerful for studying the short distance behaviour of correlators, whereas it 
is not obvious how to extract the large distance asymptotics. Also it is not straight- 
forward to extend this approach to the critical regime — 1 < A < 1 or to include an 
external magnetic field. 

These issues can be very naturally addressed in the framework of the Dual Field Ap- 
proach (DFA) to correlation functions, which was pioneered in |y], [12], [13], [Ii| ^4|, ^| 
for the example of the 5-function Bose gasHIl. A detailed and complete exhibition 
of this work can be found in the book p2[ . By means of the DFA it is possible to de- 
rive determinant representations for correlation functions in practically any integrable 
model, for which the Algebraic Bethe Ansatz can be formulated. Using the deter- 
minant representation one can then obtain explicit expressions for the large distance 
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asymptotics of correlation functions (even at finite temperature), and the inclusion of 
an external magnetic field poses no problem. The DFA thus nicely complements the 
approach of the Kyoto group. Here we will apply the DFA to the Heisenberg XXZ 
chain at zero temperature in a magnetic field h, i.e. the hamiltonian (with periodic 
boundary conditions) 

H = T,°i°s+i + "K-i + A - 1 )~ h J2 (7 j > -i< a< i, (i) 

i=i i=i 
where a z = Q a x = a y = (° and A = cos(2t/), | < 7/ < tt. 

Our discussion below follows the logic of the DFA: in section 2 we recall the 
formulation of the Algebraic Bethe Ansatz (ABA) for the XXZ chain, and then use the 
ABA to express correlation functions in terms of determinants of Fredholm integral 
operators. In section 3 we discuss how to embed these determinants in systems of 
integrable integro-difference equations (IDE) (this is analogous to describing quantum 
correlation functions by means of differential equationsBil). 



2. Determinant Representations 



Let us first review th 
XXZ Heisenberg magnetfi 
Inverse Scattering Method 
the Yang-Baxter equation^ 



i?(A,/i) 



features of the Algebraic Bethe Ansatz (ABA) for the 
la. Starting point and central object of the Quantum 



ee e.g. [Z7,|SP|) is the R- matrix, which is a solution of 



For the case of the XXZ model it is of the form 



ff(jh A) 




V o 





l 
o 



where 



sinh(A — fi+ 2irj) 




1 

o 



o 
o 



/0M)y 



i sin(2?7) 



(2) 



(3) 



sinh(A — fx) ' sinh(A — (j) 

The R-matrix is a linear operator on the tensor product of two two-dimensional linear 
spaces: R(fi) € End(<C 2 0(C 2 ). From the R-matrix (0) one can construct an L-operator 



of a "fundamental spin model" (see e.g. |22| p. 126) by considering i?(/i)n, (where n 
is the permutation matrix on (D 2 ® (D 2 ) as an operator-valued matrix by identifying 
one of the linear spaces with the two-dimensional Hilbert space 7i n of S"?7(2)-spins 
over the n'th site of a lattice of length L 



(4) 



' sinh(/i — ir\o z n ) —i sin(2?7)cr ri 
—i sin (277) cr+ sinh(yU + ir/a^ 

The Yang-Baxter equation for R implies the following relations for the L-operator 
R{\ - y) (L n (X) ® L n (fM)) = (L n (fx) ® L n (X)) R{\ - y) . (5) 
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From the ultralocal L-operator the monodromy matrix is constructed as 

Eq. ([5]) can be lifted to the level of the monodromy matrix 

R(X - n) (T(A) ® T(/i)) = (T(ai) ® T(A)) L>(A - /i) . (7) 

Below we will repeatedly use especially the following matrix elements of (0) 

[B(X),B(fi)\ = 0=[C(\),C(v)] 

iB(\),C(n)} = g{\n) (D(X)A(fx) - D(fx)A(X)) 

D(ji)B{\) = f{X,^B(X)D(ji) + g(j 1 ,X)B( Ji )D(X) 

A(fM)B(X) = f(tM,X)B(X)A{ix)+g{X,ix)B(riA(X). 

By tracing (^) over the matrix space one obtains a one-parameter family of commuting 
transfer matrices r(/i) = tr(T(ji)) = A(u) + D(u): [r(/x), t(v)] = 0. The transfer 
matrix is the generating functional of an infinite number of mutually commuting 
conserved quantum operators (via expansion in powers of spectral parameter). One 
of these operators is the hamiltonian 



d 

H = -2ism(2r]) — \n(r(n)) 

OjJL 



- 2Lcos(2r/) - 2hS z 

-it] 



Below we also make use of some properties of inhomogenous XXZ models, which are 
constructed in the following way: we first note that the intertwining relation for the 
L-operator (|5|) still holds, if we shift both spectral parameters A and \x by an arbitrary 



amount u n , i.e. 



R(X - n) (L n (X - u n ) ® L n (ji - u n )) = (L n (ji - u n ) ® L n (X - u n )) R(X - . (9) 

The reason for this is of course that the i?-matrix only depends on the difference of 
spectral parameters. We now can construct a monodromy matrix as 

U^^A-^UA-^i)---^-^)^^ ^) • (10) 

The inhomogeneous monodromy matrix ([H]) obeys the same intertwining relation (|^) 
as (|). 

The ABA deals with the construction of simultaneous eigenstates of the transfer 
matrix and the hamiltonian. Starting point is the choice of a reference state, which 
is a trivial eigenstate of t(/i). In our case we make the choice |0) = (8>£ =1 | |) n , we 
choose the completely ferromagnetic state. The action of the L-operator (H) on | |) n 
can be easily computed and implies the following actions of the matrix elements of 
the monodromy matrix 

AQi)\0) = a(/x)|0) , a(/i) = (sinh(> - in)) L , 
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D(ji)\0) = d((i)\0) , rf(/i) = (sinh(/i + z7 ? )) L , 
C(ji)\0) = , 
B(jj.)\0) ^0, 

From ( |TTD it follows that B(X) plays the role of a creation operator, i.e. one can 
construct a set of states of the form 

N 

*jv(Ai, ■ ■ ■ , Ajv) = IJ B(Xj)\0) . (11) 
i=i 

The requirement that the states (|TT|) ought to be eigenstates of the transfer matrix 
r(fj,) puts constraints on the allowed values of the parameters A n : the set {Xj} must 
be a solution of the following system of coupled algebraic equations, called Bethe 
equationM^ 

A,) ^ (12) 

d(X,) fe =1 /(Aj, Afc) i '-'-' iV - 

These equations are the basis for studying ground state, excitation spectrum and 
thermodynamics of Bethe Ansatz solvable models. For the case of the XXZ model 
with A > — 1 (the case we are interested in here) it was proved by C.N. Yang and 
CP. YamSI] that the ground state is characterized by a set of real Xj subject to the 
Bethe equations ([12]). Without an external magnetic field (h = 0) their number is 
N = L/2. In the thermodynamic limit the ground state is described by means of an 
integral equation for the density of spectral parameters p(A)l!2Ka 

2vrp(A) - / A dji K(X,fi) p{n) = D(X) , (13) 
J— A 

where the integral kernel K and the driving term D are given by 
* (* A ) = i o^w : —r^r , D(X) ~ 



sinh(/i — A + 2ir]) sinh(yU — A — 2irf) ' sinh(A — irj) sinh(A + irj) 

(14) 

Here A depends on the external magnetic field h. The physical picture of the ground 
state is that of a filled Fermi sea with boundaries ±A. The dressed energy of a particle 
in the sea is given by the solution of the integral equationHH 

e(A) - — / A da K(X, u) e(u) = 2h 2(sm(2r])) 2 

v ; 2ti J -a p v ' w vw sinh(A - irj) sinh(A + irj) V ; 

The requirement of the vanishing of the dressed energy at the Fermi boundary 
e(±A) = determines the dependence of A on h. For small h this relation can 
be found explicitly by means of a Wiener-Hopf analysis!!. For h > h c = (2 cos?]) 2 
the system is in the saturated ferromagnetic state, which corresponds to A = 0. 
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2.1. Two -Site Generalized Model 

For the evaluation of correlation functions the so-called "two-site generalized 
model" has proven a useful tool. From the mathematical point of view this is sim- 
ply the application of the co-product associated with the algebra defined by (0). The 
main idea is to divide the chain of length L into two parts and associate a monodromy 
matrix with both sub-chains, i.e. 

r( A i)=T(2 )A i)r(l )A i) ) T(i jA i) = (^ ) ) g^j) (i = l,2). (16) 

In terms of L-operators the monodromy matrices are given by 

T(2,n) = L L (ji)L L _ 1 (p,)...L n (ji) , 

T(l,/i) = L n _i(/x)L„_ 2 (/x) • • -£i(/i) • ( 17 ) 

By construction it is clear that both monodromy matrices T(z,/x) fulfill the same 
intertwining relation (Uf) as the complete monodromy matrix T(p). Similarly the 
reference state for the complete chain is decomposed into a direct product of reference 
states |0)j for the two sub-chains |0) = 1 0) 2 <8> |0)x- The resulting structure can be 
summarized as 

Ml*)\fyi = (k(j*)\0)i, A(A*)|0>i = di(^)|0>i, 

C(ji)\0)i = 0, A(^)|0)i/0, (18) 

where the eigenvalues a and d in (|ll]) are given by a(/z) = a 2 (/i)cti(/x) and d(/i) = 
d 2 ([i)di(fi). The creation operators B(fi) for the complete chain are decomposed as 
B(fi) = A 2 (fi) ®Bi(fi) +S 2 (a*) ®Di(h), which implies that eigenstates of the transfer 
matrix can be represented as 

N ni ri2 

n^oio) = EniwMH(^) 

j=l IJIjGlkGlI 

x/(Aj,A?) (B a (Af)|0> 2 ) ® (Bi(AJ)|0>i) , (19) 

where the sum is over all partitions {A 7 } U } of the set {A 3 } with card{A 7 } = rii, 
card{A 77 } = n 2 = N — n\. A similar equation holds for dual states 

N ni U2 

(oinc^) = Enn^(A;K(Af) 

j=i i,iijeikeii 

x/(A 77 ,A 7 ) (i^ldCAj)) (8 ( 2 <0|C 2 (A?)) . (20) 
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2.2. Reduction of Correlators to Scalar Products 



We are now in the position express correlation functions in the framework of the 
ABA. We will concentrate on the Ferromegnetic String Formation Probability (FSFP) 
correlation function, which is defined as follows 

m 

P(m) = (GS\ l[P 3 \GS), (21) 
i=i 

where \GS) is the antiferromagnetic ground state and Pj = |(erj + l) is the projection 
operator onto the state with spin up at site j. The physical meaning of P(m) is the 
probability of finding a ferromagnetic string (m adjacent parallel spins up) in the 
ground state \GS) of the model ([!]) for a given value of the magnetic field From a 
technical point of view this correlator turns out the simplest one in the DFAEHl. Let us 

now rewrite (pl|) in terms of the ABA. Noting that |(cr| + l) = exp — 
we find 



P(m) 



nf =1 CjXj) exp(q EZi °I°t) Uti B(X k )\0) 



(o|nf =1 c(A,)nf =1 s(A,)|o) 



(22) 



where the spectral parameters \j are subject the Bethe equations (|I2| ) for the ground 
state. This expression can be rewritten in terms of scalar products by using the two- 
site generalized model introduced above: we take the first sub-chain to contain sites 
1 to m and the second one sites m + 1 to L. We note that due to the projection 
(a = — oo) in (^) the partitions in ([TI]) and ( |2"DD must be such that there are no 



down-spins on the first sub-chain (only the partition where rii = survives), and 
thus 

1 N N N 

P(m) = — 2 (0\Y[C 2 (\ j )l[B 2 (\ k )\0) 2 n a i( A *H(A/), (23) 

a N j=l k =l 1=1 

where 

N N 

a N = (0\Y[C(\ j )l[B(\ k )\0) . (24) 

7 = 1 k=l 



2.3. Determinant Representations for Scalar Products 



In (^3|) we have reduced the computation of P(m) to the computation of a scalar 
product of the form 

N N 

S N = (0\YlC(\f)UB(\*)\0) . (25) 

3=1 k=l 

Following [^0|[n| we will show how to represent (|25|) as a determinant. Here we do 
not assume that the sets of spectral parameters {\ B } and {\ c } are the same, and we 
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also do not impose the Bethe ecmations (|T^) . From (|8|) and ( |TTD it follows that scalar 
products can be represented asli^l 

SN = T,f[^t)f[d^k)K N (\f A l , (26) 

A,Dj=l fe=l I l J ' 

where the sum is over all partitions of {X c } U {X B } into two sets {X A } and {A D }. 
The coefficients Kp? are functions of the Xj and are completely determined by the 
intertwining relation (|7J). In particular the K^s are identical for the homogeneous 
model (|5|) and the inhomogeneous model fllUp, i.e. the K N 's are independent of the 
inhomogeneities {v n } and also do not depend on the lattice length L as long as N < L. 
The reason for this is that the intertwining relations for the matrix elements A(fi), 
B(fi) , C(/x) and T>(fi) of ( |T(J| ) are the same as the ones for the matrix elements A(fi), 
B(fi), C(/i), D{ji) of(||) (see above). We will exploit this fact by considering special 



inhomogeneous models for which all terms but one in the sum in (26) vanish, and 
then represent this term as a determinant. The basic tool for representi ng scalar 
products as determinants is a Theorem due to Izergin, Coker and KorepinliSEil, which 
deals with determinant representations for the partition functions of inhomogeneous 
XXZ models constructed according to (|T0|): 

Theorem: Consider an inhomogeneous XXZ chain of even length N with inhomo- 
geneities z/j, j = 1 . . . N. Let |0) and |0) be the ferromagnetic reference states with 
all spins up and down respectively. Let i3(/i) and C(fi) be the creation/ annihilation 
operators over the reference state |0). Then the following determinant representation 
holds 

N N 

<6in*(Ai)io> = (oi nc(A,)io> 

3=1 3=1 

N N 

= {-1) N II II sinn (A Q -v k - in) sinh(A a - v k + in) 

a=l k=l 

x I [] sinh(A Q - \ p ) J] sinh(z/ z - v k ) det(M) ,(27) 

\l<a<(3<N l<k<l<N I 



where 



. . i sin(2?7) . , 

smh(A a — fk — %r l) smh(A a — v k + ir\) 



Let us now derive explicit expressions for the coefficients K^. It will be convenient 
to work with the following sets of spectral parameters 

{A AC } = {\ A } n {A c }, {A DC } = {\ D } n {A c }, 

{\ AB } = {\ A } n {\ B }, {\ DB } = {\ D } n {x B }, (29) 
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with cardinalities 

n = card{A DC } = card{A AB } , 
N-n = card{A AC } = card{A DB } . (30) 

The partition with n = is characterized by {A AC } = {A c }, {X DB } = {\ B }, 
{X AB } = = {X DC }. The corresponding coefficient K N j^cj \\ B }) * s cane d 

highest coefficient. For highest coefficients the following determinant representation 
holds 

Kn {{P} \P}) = fn^f.AfMAf.Afjjn^A^detCM^) 



h(X,ii) 



/(A,//) sinh(A — fx + 2irj) 



g(X,fJ l ) ism(2n) 

(m b ) -hx c x b ) = 9{x ^ x * ] - ^ sin ( 2 ^)) 2 m) 

{ m c] jk ^j^k) h(Xf,X B ) sinh(Af -Af + 2^)sinh(Af-Af) ' l J 

The proof of fl31~|) is as follows: Consider an inhomogeneous XXZ model on a lattice of 
length N with inhomogeneities Uj = X^ + in. We have a(A) = IIj=i sinh(A — X^ — 2in) 
and d(X) = I\f=i sinh(A — X^f). Inspection of fl26|) yields that in this situation only 
one term in the sum of the r.h.s of ( p6[ ) survives, namely the one with {X D } = {X B }. 
Thus for this special scalar product we obtain 

c • =Kn (\x c \ un)n^ h (Af-A^-2^)nsinh(A^-Af). (32) 

=A . +irj \ \ J 1 / / a i. m I 



s 



N 



On the other hand -B(A) flips one spin, and as we have chosen N to be the length of 
the lattice we find that YljLi B(Xj)\0) is proportional to the ferromagnetic state with 
all spins flipped, and thus orthogonal to all states in a basis other than |0). Thus 

N N 

Sn r = <0| J! C(Af )|0><0| n ^(Af)|0>. (33) 

V j=\<j+iri j=1 fe=1 

By the Theorem both factors can be represented as determinants. By direct compu- 
tation we find for one of the factors 



N 

<0| II C(\f)\Q) = II sinll ( A fc - Af - 2i V ) . (34) 



Using the determinant representation given by the Theorem on the other factor we 
arrive at fl3~T|). 

Arbitrary coefficients are expressed in terms of highest coefficients as follows 

W A I 1 / / \j£ACk£DC I UeABmeDB 
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/ {A AB } {A"°}\ /{A* 7 } {A Dfl }\ (35) 
\{X AB } {X DC } J n \{X AC } {X DB } J ' 

To prove fl3~5|) we consider an inhomogeneous XXZ model with inhomogeneities {uj} = 

f \X C \ { X B \ \ 

{X AB + if]} U {X AC + in}. Now only the term proportional to Kn y {A 13 } ) * n 

the sum on the r.h.s of. p6| ) survives. Proceeding as above we arrive at fl35|). 

Combining (|31|) and (|35|) with fl26|) we obtain the following expression for general 
scalar products of XXZ magnets 

Sn = n^(^>^)^(Af,Af)Esgn(Pc)sgn(P B )nMA/ B ,Ar)n^(^ C ,A^) 

j>fc j',fc l,m 

x II MAf C , Aj? °) II MAf , A^)det(M^)det(M^g) , (36) 

l,k j,m 

where P c is the permutation {A^ c , . . . , X AC , Af c , . . . , Af of {Af , . . . , A^}, P B is 
the permutation {Af B , . . . , Af B , A^ B , . . . , \#L n } of {Af , . . . , Af }, sgn(P) is the sign 
of the permutation P, and 



(M^) 7fc = t(Af,Ar)rf(Ar)a(Af), t(X,») = . (37) 



The most important step in the DFA follows next: we introduce dual quantum 
Gelds in order to simplify ( |SE| ) and obtain a manageable expression for scalar products. 
This step was first carried out for the delta-function Bose gas£9. The XXZ case can 
be treated very similarly, so that we will be brief in our discussion. The fundamental 
observation is that the r.h.s. in(p6|) looks like the determinant of the sum of two 
matrices, i.e. let A and B be two N x N matrices over (D. Then the determinant of 
their sum can be decomposed as follows 

det{A + B) = sgn(P r )sgn(P c )det(Ap r p c )det( J Bp rP J . (38) 

Pr,Pc 

Here P r and P c are partitions of the N rows and columns into two subsets 1Z, 7t 
and C C of cardinalities n (for 71, C) and N — n (for 71, C) respectively, A PrPc is the 
n x n matrix obtained from A by removing all 7?.- rows and C-columns, and B>p r p c is 
the N — n x N — n matrix obtained from B by removing all 7£-rows and C-columns. 
Finally sgn(P r ) is the parity of the permutation obtained from (1, . . . , N) by moving 
all 1Z- rows to the front. Comparison of (|36| ) with ( |3~8D shows that one does not get 



the h(X, /^-factors by simply taking the determinant of the sum of the matrices Mjk- 
This leads to the introduction of the auxiliary variables $a(A) and $p(A), which are 
called dual quantum Gelds and which are represented as sums of "momenta" Pa and 
"coordinates" Qa as follows 

MA) = Qa(X) + P d (X), $ d (X)=Q d (X)+P a (X), 
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[Pd(\),Qd(p)] = HH^til [Pa{X),Qa(»)) = HH»^)) ■ (39) 

All other commutators of P's and Q's vanish. A very important property of the fields 
$ is that they commute for different values of spectral parameters 

[$4(A), $d(h)] = = [$a(A), *a(h)) = [*d(A), • (40) 

The dual quantum fields act on a bosonic Fock space with reference states |0) and (0| 
defined by 

P o (A)|0) = 0, (0|Q a (A) = 0, a = A,D ,(0\0) = 1 . (41) 

Using the dual fields it is now possible to recast ([36|) as a determinant of the sum of 
two matrices in the following way 

Sn = Ug(\f,\°)g(\*,\f)(0\detS\0), 

j>k 

S jk = t(Af,Af)a(Af)d(Af)exp($ A (Af) + $ D (Af)) 

+f (A? , A?)d(A?)a(A?) exp (Wa?) + M*?)) . (42) 



It is possible to further simplify (^) by eliminating one dual field: we define a new 
dual vacuum (0| according to 

(6| = (0| exp [EMXf ) + P A (Xf)j , (0[0) = 1 , (43) 
and a new dual field 

<p(\) = p(A) + g(A), 

q(X) = Qa(X)-Q d (X)-(6\Q a (X)-Q d (X)\0), 
p(A) = P D (X)-P A (X), (0|g(A) = = p(A)|0), 
[p(A), ?(/•*)] = -ln(/i(A,/i)/i(/i, A)) , 

|p(A),p(A*)] = 0=[g(A),g( A i)] = ^(A) l ^(A*)]. (44) 
In terms of this field we obtain the following determinant representation 

N 

Sn = Ug(Xf,X^g(XlXf)Yl<Xf)d(Xf)Uh(Xf,X^(6\detS\0), 

j>k j=l j,k 

' jk = t(Af,Af)+t(Af,Af)^||exp(^(Af)-^(Af)) 

x -q hjXi 3 ; Xm)h{Xm> Xf) 
m=l ^(Am> Xk)h(Xf, A^) 

where r(A) = 



5 



10 



Norms of Bethe wave functions are special cases of (|25|). These were first conjec- 
tured in (see also ||). This conjecture was generalized and proved in |19[. The 
result is 



N 

(oi n aw n B{x k )\o) = n x k ) n a^MA^det w 

j=l k=l j^k j=l 



N 

n 

fc=i 



N 



M'jk = sm(2r]) -K(Xj, X k ) + i S jk 



sin(2r/) [-K(Xj,X k ) +5 



_d_ 

dX j [ 
d 



dX, 



ln(r(A,)) + J2 Hjf^) 

N 1\ 

■Hr(Xj)) + Y,K(X 3 ,Xn) , (46) 

n=l J / 



where K(X,fi) and h(X, fi) are defined in (PI). We note that (|46| ) can be obtained 
directly from fl45|) by setting all dual fields to zero and then taking the sets of spectral 
parameters {X c } and {X B } equal and imposing the Bethe equations (|12|). 

2.4- Determinant Representation for FSFP correlation function 

Let us now use the machinery built up above express P(m) as a determinant. We 
will proceed in two steps: we first will analyse ( [23D without using that {X B } = {X c } = 
{A} and without imposing the Bethe-equations (|T^) . In the second step we will then 



impose these two constraints. Using ( [45]) we can represent the scalar product in the 
two-site generalized models in (|23|) as a determinant and obtain 



P(m) = — II^f,A^(Af,Af)(0|det S ({A c },{A B })|0) 



a N 



(s({x c },{x B m k = ai(AyH(A£) 



j>k 
C\ 



t(Xf, X B )a 2 (Xf)d 2 (X B ) exp ($ A (Af ) + $d(A 



+t(Af , Af )d 2 (Xf )a 2 (X B ) exp ($ D (Af ) + $a(A 



(47) 



Here the dual fields are defined according to 

MA) = Qa(X) + P d (X), $ d (A)=Q d (A)+P a (A), 

[p d (X),q d (v)} = HKKri), [P A (x),Q A (ri} = HKv,x)) . 



(4* 



All other commutators vanish. The reference state |0) and its dual (0| are defined via 

P a (A)|0) = 0, (0|Q a (A) = , a = A,D ,(0|0) = 1 . (49) 

So far we have not used the fact that we are dealing with expectation values of 
Bethe states, i.e. we have neither used the fact that {X c } = {X B } = {A} nor imposed 
the Bethe equations fljjl). Let us now impose these constraints. Like for the case of 
scalar products one of the dual fields can be eliminated 
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We define a new dual vacuum (0| and a new dual field according to 

(0| = (0|exp rtp D {\j) +P A (X j ) ) j , (0|0) = 1 , 

0(A) = $ A (A) - $ D (A) , (50) 

The momenta p (A) and coordinates q(X) of the dual field 0(A) obey the commutation 
relations 

[q(li),p(\)]=hi(h(\,fi)h(fi,\)) (51) 

(all other commutators vanish). By straighforward rewriting of (|47| ) in terms of the 
new field and the new dual reference state we obtain 

1 N 
P{m) = — n/(A J) A fc )/(A fc) A,)n«(A i )rf(A,)(0|detA ; i|0) , 

a N j>k j=1 

M jk = t(A J ,A fc )+t(A fe ,A :; )^4exp(0(A fc )-0(A J )). (52) 

r 2{Aj) 

Here we have used that 

(0|exp(f:<|. A (A J ) + $ D (A J )] =nMA i; A fc )(0| . (53) 

It is found that whereas p(A)|0) = 0, the coordinate g(A) of 0(A) do not annihilate 
the new dual reference state (0|. Therefore we "shift" 0(A) by subtracting its vacuum 
expectation value in analogy with fl44j) 

<p(\) = 0(A) - (O|0(A)|O) = p(A) + g(A) . (54) 



By construction p and q have the same commutation relations (|54|) as the mo- 
menta/coordinates p(X) and g(A) of 0(A). Furthermore p(A)|0) = and (0|g(A) = 0. 
The shift is equal to 

K( A) = (0 W A)|0)=p°(^) . ,55) 

If we replace in (|5*2"| ) by ip we pick up additional factors due to the shifts 

1 N 
P(m) = — n/^A^n^KAjOCOldetCflO) , 

°n j+k j=1 

Qjk = t(A J ,A fc )+t(A fc ,A i )^^e^)-^)e K ( A ^- K ^). (56) 
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The off-diagonal matrix elements of G can be further simplified by simply imposing 
the Bethe equations. Rewriting the Bethe equations ( |T2| ) as 



j=i X k ) 

i N 



n(x h 



n 



h(Xj, \i 



\N-1 



r 2(A fc ) ~ =1 h(X k , Xj) 



T\{Xh 



N 



(57) 



we find that the additional factors take the form 

r 2 (X k ) _ ri(Aj) 



exp(K(X k ) - k(Xj)) 



rziXj) n(x k ) 



(5? 



To get the diagonal matrix elements we have to investigate the limit Xj — > X k of 
(|56|) in detail. In the limit Xj — > A& the sum of the first two terms in Gj k and the 
expression in brackets are both of the form "~" . By using l'Hospital's rule we find 



lim 



t(Xj, Afc) + t(X k , Xj) 



+ sinh(2^) 



rjjjXk) 
r2(Xj) 



OX 



d 

+ sinh(2i?7) — 

A=A, OA, 



Hr 2 {Xj))+Y,H 



—2 cosh(2z7/) 
h(Xj, X n 



h(X n , Xj) 



• (59) 



Putting now everything together we obtain the following representation of the FSFP 
as a ratio of two determinants 



P(m) 
Gj k 



(6[dgtg|0) 

detJV' ' 
t(Xj, X k ) + t(X k , Xj 



r 1 (A j 



exp ((f(X k ) - <p(Xj)) 



n(X k 

+5 jk sm(2 V ) \LD(Xj) + £ K(X j7 X n )) 



(60) 



where (p and A/ 7 are defined in ([54]) and (^) respectively and where D and K are 
defined in (0). Eqn ( |60| ) is our final result for the FSFP on a finite chain of length 
L. As always in the Bethe Ansatz things simplify essentially in the thermodynamic 
limit. Let us first investigate the thermodynamic limit for the norm (^6[). We 
write A/ 7 as the product of two matrices: 



N'j k = sin(277) tjmJi 



K 



mki Ijm 3jm n j Jjm &jrrfirm (^"0 
m "m 

where 9 m = LD(X m ) + J2 n K(^m, X n ). The determinant of A/ 7 is the product of the 
determinants of / and J. Next we use that the set of roots {A.,} describes the ground 
state and the roots thus obey the equations 

27rLp(Xj) - J2 K (*j, **) = LD(Xj) , j = 1 ... N , (62) 
fe=i 
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which is the discrete version of (|13|). Here p(Xj) = L ^ A . \_ X ) i which becomes p(A) 
defined by ( |13"D in the thermodynamic limit. We thus can rewrite 9 m = 2nLp(X m ), 
which leads to 

N 

det J = [] 27rLp(A j ) . (63) 

3=1 

In the thermodynamic limit the matrix I turns into an integral operator I = id—^K 

1 f A 



I*f 



/(A)-—/ dpK(X,p)f(p) , (64) 

271 J-A 



where K is the kernel of K defined by (14). The matrix Gj k in ( |60D is treated in a 
very similar way. We rewrite it product 

G jk = sin(2r]) W jm J mk , (65) 

m 

where Jj m = 6j m 27iLp(X m ) is the same as above, and 

W jk = S 3 k + . * U\j, X k ) + t(X k , A,)^4 exp (^(A fc ) - ^(A,)) ) . (66) 

In the thermodynamic limit the matrix Wjk turns into an integral operator W = 
id + ^-V with kernel 

Z7T 



V™(\p) = ~ Sin(2r/) ( - e-\X)e(p) | 

sinh(A — /i) lsinh(A — p + 2irj) sinh(/i — A + 2ii]) J 



\sinh(A — 

Thus in the thermodynamic limit the FSFP P(m) for the XXZ magnet can be rep- 
resented as a ratio of determinants of Fredholm integral operators (id + -^V^j and 

(id — -^Kj in the following way 

= (0|dt(tt + ^)|0) 
det (id -i-ST) 

Here (0| and |0) are the vacua of the dual bosonic Fock space defined in ([501). 



3. Integro-Difference Equations 

In this section we embed the numerator of the r.h.s. of fl6"9|) in a system of integro- 
difference equations. The denominator is independent of the distance m and thus 
merely amounts to an overall normalization of P(m), which is difficult to determine 
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anyhow and will be dropped in what follows. We start by bringing the kernel of 
V ([J7|) to "standard" formE3. This is done in several steps. We first perform the 
similarity transformation S with kernel S(X) = (e(A))a on 1 + ^V, which leaves the 
determinant unchanged. This "symmetrizes" the kernel 

V {m) (X a) = ~ sin ( 2r ?) ( e-\\)e(fi) \ 

' sinh(A — jj) lsinh(A — fi + 2ir]) sinh(/i — A + 2iif) J ' 

where 

/sinh(A+^)\f M X \ cti\ 
e(A)= Uinh(A-^) J ex P^)- ( 71 ) 

Next we change variables to /3 — exp(2A), a = exp(2/i), w = exp(2ir]). Elementary 
calculations yield 

pw - l\ 2 /aw - 1 
p-w ) \a-w 

x ij?m-«m- a ^x (72) 



Next we use the identity 

roo / a \ — i 

/ ds exp -is(pw ) = , (73) 

Jo V w J pw — — 

(which holds as a > 0, P > and sin(2?7) < 0) in order to eliminate the unwanted 
denominators in (|72|), and finally perform yet another change of variables to z(ct) = 
. This change of variables maps the real axis on the contour C : a — > z = exp ia 
where — ip <a<27r + ?/><?7 (■?/.> <0 by definition). The endpoints £ = e 1 ^ and 
£ = e~ 1 ^ (we integrate from £ to £) of the contour are related to the magnetic field h 
and the anisotropy A. After elementary computations we obtain 



da^^ln(P) ln(a). dz 2 fwzi — 1 Z\ — w\ r°° ^ e+\zi\s)e_(z 2 \s) — Z\ <-> z 2 



2a 2 2 i \wz 2 — 1 z 2 — w) Jo Z\ — z 2 

where the functions are given by 

ef\z\s) = (2sm(2r ] )^^z- m e-^A\x P (i- Sir: ~ 



(74) 



z — w / \ w z — w 



(75) 



e ( : n) (z\s) = (2sm{2ri) ^—^z m e^ z) Yexp(-tsw— 



W J \ z — w 
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Note that the "momenta" p(z) and "coordiantes" q(z) of the dual field <fi(z) now obey 
the commutation relations 



[q(z 2 ),p(zi)] = In 



1 (WZ 2 — 1 Z\ — W WZ\ — \z 2 — W 2 _ 2 

H w — w 



4 sin 2 (2r]) Vw^i — 1 z 2 — w ' W2 2 — 1 — u> 

(76) 

Eliminating the factor in brackets before the integral in ( |74| ) by a similarity transfor- 
mation we obtain 

det(l + ^-V) = det(l + y (m) ) 

(w) / I \ (wi) / i \ (m) / i \ (m) i i \ 

v M {zilz2) = J_ f ds elWH W)~eL WKW), 
2n Jo zi — z 2 

The integral operator now acts on a function f(z) as 

(V^f) ( Zl ) = I dz 2 V^( Zl \z 2 )f(z 2 ), (78) 

where the integration is to be performed along the contour C. We note that V is 
symmetric and nonsingular at Z\ = z 2 . Eqn ( |T?D is the desired standard form for the 
kernel, and was obtained for the price of introducing the auxiliary s-integration. 
The resolvent of is defined by (1 + V^)(l - R {m) L= 1 and its kernel 
(zi\z 2 ) can be written in a form similar to Eq. (|TTD, namely 

= — i. r is i^i^-i^MftHM*) (79) 

2it Jo Zi — z 2 

Here the functions f± are solutions of the linear integral equations 

f ± ( Zl \s) + [ dz 2 V( Zl \z 2 )f ± (z 2 \s) = e±( Zl \s) . (80) 

In terms of these functions we introduce the integral operators B^\ a, b — ± acting 
as (B^f^)(s) = J °° dt B%\s, t)f{t) with the kernel 

B%\s,t) = ^f c ^fjr\z\s)et n) (z\t), a,b = ±. (81) 

The transpose B T acts like [b j fj (s) = / °° dt B(t, s)f(t). We are now able to 
formulate the embedding of P{m) into a set of integrable integro-difference equations: 

(i) The lattice logarithmic derivative of det (l + vw) is given in terms of the 

det (l + V^ m+] l ' 



integral operator B^ as 



det (l + V( 



III) 



det (l + B {m ^) (82) 
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(ii) The logarithmic derivative of det (l + V^ m ^ with respect to the boundaries of 
the contour C is expressed in terms of the functions F± (s) = /£ (£|s) and 
G ( r\s) =/i m) (e I s) as follows 

-z^lndet(l + t/) = i- j™ ds {F + (s)d^(s) - F_{s)d^F + {s) 

+G.{s)d^G + {s) - G+(s)fyG_(s)} 

QHds F + ( S )G_( S ) - F_( S )G + ( S ))83) 

(hi) The following set of completely integrable integro-difference equations for the 
unknowns F, G, B in (i) and (ii) holds 



1 P (m+1) 
1 P (m+1) 



' ' — i(m) 



l:T ] - ( l B 



(m+l)\ ( B {m)\ T p (m) 



(m+l) f D (m)A T> \ n(m) _ ^(m+1) 



j=G? +1) = i{GV m) -(l-5t +1 0K+) TG(m) } 



1 + ^)^4^ 

(84) 



1 ^(m+l) 



1 | + B (™+D ( 5 (_™)) T ) G w - B {m _ +1) (l + (5fJ) T ) GV m) } 



and 



i^rB ab { S ,t) = ^-lF a (s)iF b (t) + (F + B_ b -F_B +b )(t)} 



+G a (s) [G b (t) + (G+B_ fe - G_B +{) ) (t)] | (85) 



where a, b = ±. 

The proof of (i)-(iii) is analogous to the one for the XXX-case so that we only sketch 
the main steps, (i) is a direct consequence of the shift- equation 

2tt Jo 



Z\ 



>6) 



which follows directly from (|67D and (76), whereas (ii) and (hi) are consequences of 
the Lax-representation 



^f { + m+l \z\s) = -fT\z\s)- 
y z z 



(m). 



1 - B 



(m+l) 



37) 
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and 



-^fMs) = ^\^^m\s)fdt(umfMt)-f-mu(m) 

+j^f±mfdt (umf-(z\t) - urn u^\t))\ m 

This concludes our derivation of the embedding of P(m) in a system of integrable 
(due to the Lax-representation (|87| ) - (|88|) ) integro-difference equations. We note that 
there is a integral-operator-valued Riemann-Hilbert problem associated with these 
IDeI, and that the asymptotics of P{m) for large m can be extracted from the 
solution of the Riemann-Hilbert problem. So far this has been done only for some 
limiting casesi, but we hope to report the full solution in the future. 
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